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Abstract 

We describe a new class of spin liquids with global SU(2) spin rotation symmetry in spin 1/2 
systems on the triangular lattice, which have real Majorana fermion excitations carrying spin S =1. 
The simplest translationally-invariant mean-field state on the triangular lattice breaks time-reversal 
symmetry and is stable to fluctuations. It generically possesses gapless excitations along 3 Fermi 
lines in the Brillouin zone. These intersect at a single point where the excitations scale with a 
dynamic exponent 2; = 3. An external magnetic field has no orbital coupling to the SU(2) spin 
rotation-invariant fermion bilinears that can give rise to a transverse thermal conductivity, thus 
leading to the absence of a thermal Hall effect. The Zeeman coupling is found to gap out two- 
thirds of the z = 3 excitations near the intersection point and this leads to a suppression of the low 
temperature specific heat, the spin susceptibility and the Wilson ratio. We also compute physical 
properties in the presence of weak disorder and discuss possible connections to recent experiments 
on organic insulators. 



I. INTRODUCTION 



The recent experimental evidence for spin liquids in the triangular lattice organic com- 
pounds K-(ET)2Cu2(CN)3p^and EtMe3Sb[Pd(dmit)2]J^has sparked much interest in char- 
acterizing the experimental signatures of the many candidate spin liquid states. 

For the compound fi;-(ET)2Cu2(CN)3, a theorj^ building upon the proximity of a magnetic 
ordering quantum critical point is compatible with the recent observation of magnetic order 
induced by a small external fielcP. 

On the other hand, EtMe3Sb[Pd(dmit)2]2 is characterizecP by a thermal conductivity, 

for which k/T reaches a non-zero limit as the temperature T — )■ 0, and this is strong 
evidence for the presence of gapless excitations across a Fermi surface. A spin liquid state 
with a spinon Fermi surface has been proposecP^, and so is a natural candidate for this 
material. However, this spinon Fermi surface state is also expected to display a thermal Hall 
effect!^ and this effect has not been detected so fai^. 

This paper will examine another possibility for a spin liquid state with a Fermi surface 
of spin 1 excitations. We will assume that the Fermi surface excitations are real Majorana 
fermions and this, as we will see, allows us to retain the longitudinal thermal conductivity 
while suppressing the thermal Hall effect. 

Our approach relies upon following representation of S* = 1/2 spins in terms of = 1 
Majorana fermiond^^HE] 

S^" = le^"''7V- (1-1) 

Here we have suppressed site indices, and the Majorana fermion operators all anti-commute 
with each other, and have a unit square (7")^ = 1 (no sum over a). As explained by 
Shastry and SerP^, such Majorana fermions provide a redundant but faithful realization of 
the Hilbert space of S" = 1/2 fermions. The redundancy is linked to a Z2 gauge invariance 
7" — )■ —7", which then also plays a crucial role in the description of any spin liquid states; 
some related issues are discussed in Appendix |X] 



The representation in equation (1.1 ) has been used extensively in recent workP^, follow- 
ing the exactly solvable spin model proposed by Kitaev'^. A rich variety of solvable models 
have been found on different types of lattices, some with global SU(2) symmetrjJ^^ESI^ oth- 
ers with Fermi surface^^MIESI However, none of them are on the triangular lattice, and 
none of them have both SU(2) symmetry and a Fermi surface: these are clearly important 
requirements for making contact with the experiments on EtMe3Sb[Pd(dmit)2]2- 

Here, we shall not attempt to find an exact solution to a particular model Hamiltonian. 
Instead, we will build upon the extensive experience that has been gained by parton con- 
structions of mean-field spin liquid states, and the establishment of their stability by an 
effective gauge theory of fluctuations. Previous constructions of Z2 spin liquids relied upon 
writing the spins either in terms of Schwinger bosonP^ or fermion^, and here we will apply 
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an analogous analysis to the Majorana parton construction in equation (1.1). We will be 
aided in this analysis by the Projective Symmetry Group (PSG)'^ which we shall apply to 
the effective Hamiltonian for the Majorana excitations. 



A. Low energy theory 

We begin by postulating the existence of a SU(2) invariant spin liquid state on the 
triangular lattice, whose quasiparticles are described by a triplet Majorana field 7°(t'), 



a = x,y or z. Although we are using the same notation as in equation (1.1), the Majorana 



field operators used in the low energy field theory create the physical quasiparticles and 



so can be strongly renormalized from the underlying Majorana fermion in equation (1.1). 
Noting that the Majorana bilinear Hamiltonian has to change sign both under time reversal 
(TR) and under a lattice rotation by vr, we assume that the 7°" transform trivially, i.e., 
7" — )• 7" (without a possible sign change) under all the PSG operations associated with a 
modified triangular lattice space group. In this modified triangular lattice space group the 
elementary operation of rotation by 7r/3 is replaced by the same operation compounded with 
TR. Furthermore, the Majorana operators transform naturally in the S = 1 representation 
of spin rotations and are real operators which are invariant under time reversal. These simple 
and general transformation rules are already sufficient to strongly constrain the effective low 
energy theory of the 7"(r). 

Let us begin by writing an effective Hamiltonian for the 7° bilinears as an expansion in 
spatial gradients. 

Demanding invariance under by 27r/3 rotations (a double application of TR + tt/3 rota- 
tion) and hermicity, we are led to the Hamiltonian 

Ho = two J d^r 7" (Pi +V2 + P3) T (1.2) 

where the Vi = 5i ■ 'V are directional derivatives along the 3 principal directions <5i, 82 and 
83 shown in Figure |2] and wq is a real number. We remark here that there is no term without 
spatial gradients because the Majorana fermions square to unity and because of SU{2) spin 
rotation symmetry. However, we clearly have 61 + 62 + 63 = and so Hq vanishes identically. 
To obtain a non-zero contribution, we have to expand all the way to 3 derivatives'^ , when 
we obtain two independent terms which can be written as 

n = i j [w,V,V2V3 - W2 {VIV2 + VIV3 + VlV,)] 7" (1.3) 



where both wi and are real parameters. The low energy Hamiltonian in equation (1.3) 



underlies all the results derived in this paper. From this it follows that the long wavelength 
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FIG. 1. The Fermi surface (shown as black hnes) corresponding to the lattice dispersion in equa- 
tion (2.5) in Section |ll] — equation (1.3) is the generic continuum version of the same theory. The 



BZ is bounded by the green border while the 'occupied' states are shaded. 



excitations of this theory have the dispersion 



Eq'^^\q^cos{39q + ^) (1.4) 
t cos = |(2wi + W2) 



where 



tsin0 = 



Next, let us describe the structure of the low energy excitations oi l-L. As we shall 
demonstrate shortly in Section [TT] and as is illustrated in Figure [T| there are two classes 
of excitations. First, there are the excitations with momentum qr 0, which have energy 
~ \q\^i and so look like those of a quantum-critical theory with dynamic exponent z = ?>. 
Second, there are the linearly dispersing gapless excitations along Fermi lines which meet at 
g = 0. 

It is now straightforward to establish the perturbative stability of "H. The collective modes 
arising from decoupling spin interaction terms like the four Majorana exchange interactions 
in equation (2.1) constitute the gauge fluctuations of a Z2 gauge theorj^^^Msi xhese have a 



finite range of stability without a transition to confinement, when our theory as described 
above is valid. Next, we consider the influence of terms quartic in the Majorana fermions. 
These quartic couplings will lead to innocuous Fermi liquid renormalizations of quasiparticles 
along the Fermi lines, just as in any Fermi liquid. The influence of quartic couplings on the 
2; = 3 excitations near qr = is more subtle, but can be analyzed by a standard scaling 
argument. The scaling dimension of 7° is d/2, where c? = 2 is the spatial dimension"^ . 
Dimensional analysis now shows that a quartic coupling with p spatial derivatives has scaling 
dimension z — d — p. For z = 3 and d = 2, this is irrelevant only if p > 1. With the 
requirements of SU(2) invariance, it is easy to show that any quartic term must have at least 
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p = 2 derivatives: the simplest non-vanishing term with SU(2) spin rotation invariance has 
the generic structure ~ (7"c}7°)(7'^97'^). We emphasize that SU(2) spin rotation symmetry 
is crucial to the stability of the theory: in its absence, marginal quartic terms arise that 
could destabilize the postulated liquid. 

The outline of the rest of this paper is as follows. In Section [TT] we shall derive the 
Majorana mean field theory corresponding to the Heisenberg antiferromagnet and reproduce 



the general low energy spectrum postulated in equation (1.3). We shall also demonstrate 



how any Majorana bilinear Hamiltonian on the triangular lattice gives rise to equation (1.3) 



in the low energy long wavelength limit. We will then describe the experimentally observable 



properties of this state in Section |III[ The influence of weak disorder will be presented in 
Section |IVl 



II. THE MEAN FIELD MAJORANA HAMILTONIAN ON A TRIANGULAR 
LATTICE 

A. Majorana mean field theory from a spin Hamiltonian: an example 

In this section we derive a Majorana mean field theory starting from the AF Heisenberg 
model on a triangular lattice, mainly to demonstrate the mechanics of such a derivation^ . 
We choose a real space coordinate system such that one set of bonds point along the x-axis, 
as shown in Figure [2] The Hamiltonian has the form: 



J 



J2 l''{x)-f''{x + S)-f^{x)-f^{x + S) (2.1) 



where 5 are the three nearest neighbor bonds labeled in Figure [2] which are related by 
rotations through 27r/3. We perform a mean field analysis for the above Hamiltonian by 
assuming the mean field ansatz'^ 

{nxh^x + 6))=tg5^p (2.2) 
which is also graphically represented by directed bonds in Figure [2] Using this ansatz, 



equation (2.1) becomes: 



^MF = ^ E [2 X 2zgrix)rix + 6) + 2g^] 



x,S,a 



^ ^ (7"(^)7"(^ + S)- r{x + S)r{x)) + ^-AfJg' (2.3) 

x,S,a 
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FIG. 2. The triangular lattice showing the lattice vectors and the 5 vectors used in the text. The 
arrows on the bonds define the mean field ansatz equation (2.2) that also becomes the scheme 
for assigning the same phase to hopping parameters. The amplitudes of hopping processes in the 
directions opposite to the specified bond directions pick up an additional factor of -1. 



where M is the number of sites, assumed to be an even number. This Hamiltonian can 
be diagonalized using the momentum states defined over the Brillouin zone (BZ) of the 
triangular lattice: 



= 7^ E 7'^(^)e"^''-^ ^ = V E ^^'"'^ (2-4) 

^ X q 

These h operators are complex Fermions with 6_q = 6^ and {6^,6^} = SapSp-q^ showing 



that independent h operators cover only half the BZ. Using these in equation (2.3), we can 
diagonalize the Hamiltonian: 



MF — 

a 



J2 [-JgJ^'inq-d] bl^b'^^ +^-MJg' 

.qeBZ \ q J 

2^ ^^-A +4-^^^ (^q = 8J^smysm sm 1 

J2 +lAfJg' {using E^q = -Eq) 

BZ'\E>0 ^ 



a,qeBZ'\E>0 



a,qeBZ 



Here we have introduced the notation BZ' to denote that half of the BZ where the quasipar- 
ticle energies of the spin rotation-invariant Hamiltonian are positive. The fermion creation 
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operators in BZ' are related to those in the remaining half of the BZ by the particle-hole 
relation 6_q = 6^. The structure of the Fermi sea obtained above is shown in Figure l| 
where BZ' consists of the un-shaded regions of the BZ. It follows that near qr = 0, the 



quasiparticle energy has the same form as derived earlier in equation (1.4) using a gradient 
expansion 

9^0 f Jg\ 3 



Eq^^ j q-" cosWq (2.6) 
The ground state energy ^^^g (fi' "~ f ) is minimized when g = ^"^ : 

Eo = -^JAf = -0.22 JAf (2.7) 

As expected, Eq is higher than the numerically calculated ground state energy of about 
—0.54 J per sit^^for the best candidate spin-ordered ground stat^^. Additional interactions 



should be added to the Heisenberg Hamiltonian in equation (2.1 ) to stabilize the spin liquid 
state described by the non-interacting Majorana ground state. 



B. The general low energy effective theory on the lattice 



This subsection will give an alternative presentation of the ideas of Section lA, working 
directly with the lattice Hamiltonian, rather than the continuum theory. 

A general spin SU (2) rotation-invariant and translation-invariant low energy effective 
Hamiltonian of Majorana bilinears has the form: 

-Hmf = I J] t(d)7"(a;)7"(a; + d) (2.8) 

where Hermicity requires that t{d) is real and antisymmetric in the hopping vector d: 

t{-d) = -t{d), t{d) E M (2.9) 

If this Hamiltonian describes a spin liquid, the observable quantities which are the spin 
correlation functions should not break any lattice symmetry, in addition to the spin rotation 
and lattice translation symmetries discussed above. However, since a lattice rotation by tt 



and time reversal separately flip the sign of the mean field Hamiltonian equation (2.8), the 
theory can be invariant only under a combined application of the two. This uses the fact 
that the Majorana operators are hermitian and also that due to spin rotation symmetry, 
their bilinears cannot acquire any additional factor under symmetry operations. We require 
the maximum possible adherence to the lattice point group symmetry consistent with these 
observations — a lattice rotation by 7r/3 combined with time reversal must leave the Hamil- 
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tonian invariant. This, along with invariance under lattice translations and reflection about 
a bond, are the elementary symmetry operations that deflne the class of effective Hamiltoni- 
ans which may possess the Majorana spin liquid ground state described in this work. Even 
with this reduced set of symmetry operations, all equal time correlation functions with an 
even number of spin operators will remain invariant under the full set of lattice symmetry 
operations. 



These arguments motivate rewriting the Hamiltonian in a manner that makes the anti- 
symmetry under a rotation by vr/S apparent: 

5 

%MF = t i(^i)$^(-l)"7"N7"(^ + ^W3(d)) (2.10) 

x,{d},a n=0 

where {d} denotes the set of hopping vectors, modulo those that are related through ro- 



tations by multiples of vr/S. In terms of the momentum state operators equation (2.4), we 
have: 

nMF = iY. (2-11) 



where 

£^rf(g) = 16t(d)sm^— sm sm (2.12) 

is the contribution to the Majorana dispersion from the hopping processes characterized by 
the hopping vector d. This expression tells us that the dispersion of the long wavelength 
low energy modes near qr = have the form 

Eq t cos(30q + 0) (2.13) 



where t and (p are real constants. These parameters are the analog of those in equation (1.4) 
obtained from corresponding continuum analysis. Figure [3] shows the Fermi sea and Fermi 
surface corresponding to a model with a next nearest neighbor hopping amplitude that is 
one-flfth of the nearest neighbor hopping amplitude and demonstrates the existence of three 
Fermi curves intersecting at g = 0, in this system. 



From this calculation it is clear that the low energy mean fleld theory is composed of two 



kinds of excitations which are smoothly connected to each other, as was noted in Section I A 
Excitations near q = have a dispersion that varies as and a dynamical exponent z = I: 
while those along the Fermi curves behave like the excitations of a 2D Fermi gas. 



8 




FIG. 3. The Fermi surface (black curves) corresponding to equation (2.10) with a next-nearest- 
neighbor hopping amphtude that is 20% the nearest neighbor hopping amphtude. The BZ is 
bounded by the green border while the 'occupied' states are shaded. 



C. The low energy effective theory in the presence of a perpendicular magnetic 
field 



Let us first consider orbital coupling terms which do not violate the SU{2) spin rotation 
symmetry. In this case the arguments that lead to the formulation of the Hamiltonian 
will be no different that in the field free case considered in the previous sections and so the 
Hamiltonian will be invariant under rotations by 27r/3. The transverse thermal conductivity, 
which involves averaging the product v^Vy over momentum space, will be zero since the sum 
of VxVy over points related by 27r/3 rotations is zero. Equivalently, it may be noted that 
the PSG implies that there is no orbital coupling between the applied magnetic field and 
a fermion bilinear: it is not possible to find a fermion bilinear which is invariant under 
translations and by spatial rotations under 7r/3. It follows that orbital coupling of the 
magnetic field will not induce a thermal Hall effectP' in our theory. This is in contrast to 
what happens in the case of the U(l) spin liquid with a spinon Fermi surfac^, where the B 
field does couple to fermion bilineaiOS the coupling is of the form S ■ (V x J), where J is 
the U(l) spinon current. 

The other way in which a perpendicular magnetic field Bz enters the Hamiltonian is via 
terms that break the spin SU{2) symmetry down to a U{1) symmetry of rotations about 
the direction of the magnetic field. Such a coupling will not affect the 7^ fermions but will 
couple the 7^'^ Majoranas into 5"^ = ±1 excitations. The most relevant term in that case 
is the Zeeman term —iB'^^'^y/2 which does not break the three-fold rotation symmetry and 
thus does not lead to a thermal Hall effect. 
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D. The spectrum in the presence of the Zeeman coupUng 



The Zeeman term —/j^qBSz does not affect the spectrum of the 7^ fermions because they 
carry spin Sz = 0. The Hamihonian of the 7^'^ fermions, however, is modified: 

= 1 Yl {Ek-sfioB)cl{k)cs{k) 

k,s=±l 

= ^{Ek + fioB)cl{k)c^{k) + c-number (2.14) 
fc 

where the new fermionic excitations with spins Sz = s/2, s = ±1 are 

c.{k) = ^I±j^; c,. = c^_, (2.15) 

The Fermi surface now consists of arcs in three of the six wedges partitioning the BZ, as 
shown in Figure |4} 




FIG. 4. The Fermi sea of c_ fermions (shaded regions) in the presence of a Zeeman term (assuming 
HqB > 0). The hexagon is the BZ for the triangular lattice. The 7^ fermions are not affected by 
a magnetic field in the z direction and will retain their original excitation structure as shown in 
Figure [l} 



10 



III. PROPERTIES OF THE CLEAN MAJORANA SPIN LIQUID 



The bilinear Majorana Hamiltonian which will be used in the following sections is 



nMF = iY.^^%K (3-1) 



2 



where the g — i- form of Eq is given by equation (2.13). The propagatoi'^ for the Majorana 
excitations is given by'^ 

W> = ^%^=^"/^^.--.^.W (3-2) 

z — niq 



A. The low energy density of states (DOS) 

Because of the dispersion, the contribution to the density of states from the states near 
A; = diverges as the energy — )■ 0. The divergence may be calculated from the effective 



Hamiltonian in equation (2.13): 



p{E) = 3 ^ 5(E - tk^ cos 3^fc) 

k 

Jo 47r^ J_^/6 

3 dkk 1 3 dkk 



2|i|7r2 7|s/t|i/3 k^ Ji _ 2\t\'K^ J\E/t\^ri ^k^ - (E/t) 



^0.18(t2|E|)-i/3 = Po|^r'/' (3.3) 
where A ~ 1 is the upper cutoff for the momentum integral. 



B. Specific Heat 

The specific heat, as T — is given by 
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C. Magnetic susceptibility 



Only the SU (2) spin rotation symmetry-breaking Zeeman term — |yUo-B"7'^7^ will give rise 
to a magnetic moment due to the application of a magnetic field Bz. The static susceptibility 
may be calculated from the spin correlation function: 



47r2 ^ {iuj^ - EkY 

lUJn 

'■^«?M2^~<i££-./3(_^) ^0.38^„p„r-'/3 (3.5) 
D. The Wilson ratio — comparison with a 2DEG 

For a spin 1/2 free fermion gas, the susceptibility and specific heat are given by 

£'2DEG TT^ 
Xxx = ^ , Cv = y£»2DEG-r (3.6) 

Thus, the Wilson ratio of this model is 

0.38poT-V3 fg^^^oT ^ 

1.18poTV3 « -^-^ ^'^•^^ 

times that of the free spin 1/2 electron gas, assuming fiQ = fiB, the Bohr magneton. 

A spin-1 non-interacting 2DEG has a Wilson ratio that is 8/3 = 2.67 times that of the 
spin-1/2 free 2DEG, i.e, smaller than that of our model. 



E. Static Structure Factor 



The spin static structure factor is given by 



(S^MS^m = --^e'^"V"'^(7''(r)7^(r)7'^(0)7'^(0)) 
16 



foru^i/ 

(3. 



(7"(r)7"(0))^ 
4 



for ijl — u (no sum over a) 



This simplification occurs due to the absence of correlation between Majorana fermions 
of different flavors arising out of spin rotation invariance; (7"(t")7'^(0)) = f{r)6ai3, where 
the function /(r) may be found as follows. Expressing the position vector r — riiai + ^202 
in terms of the lattice displacement vectors ai_2 as well as the wave vector r — piKi -\-p2K2 
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FIG. 5. Variation of (7r/(r))2 = Air^ (5"(r)S"'(0)) on the real space lattice, at T = with only 
nearest neighbor hopping. The red circle at the center is r = and there is no sum over a. We 
see that (S°(r)5"(0)) decays as an inverse square of the distance (the strongly directed nature is 
an artifact of nearest neighbor hopping) . 



in terms of the reciprocal lattice vectors Ki^2 defined through Ki ■ aj = Sij, we find that 
fif) = (7"(^)7"(0)) (iio sum over a) 



27r2 

ni,n2,ni—n2y^0 



P,Q P 
T=0 1 j j gi(pini+p2n2) 

J Jbz' 

„-i7r(ni+2n2) [ 

^ (_1 + e*-"i) (1 + e-K+2n2)) ^2 

/7r^nin2 {—ni + n2) [ 

+ (1 - e^™2) (1 + e*"("i+2"2) (1 - e^""^ + e^""^)) rii 



(3.9) 



From the last expression, which is valid only for the case with nearest-neighbor hopping, we 
can prove that at T = f{r) is zero when ni, n2, ni — n2 7^ 0. Thus, the function is non-zero 
only along 3 lines defined by ni, n2, ni — n2 = 0. The spatial variation of the squared value of 
this function at T = 0, which is proportional to the static spin correlation function, is shown 
in the Figure |5| We find that static spin correlations are negative and decay according to 
the inverse square law along the six directions (3 lines) discussed above. Inclusion of longer 
range hopping processes will modify the highly directional nature of the correlations. 
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F. Effect of a perpendicular magnetic field 



In the presence of a Zeeman term —n^BSz which is small in comparison to the bandwidth, 
the DOS gets modified to 



Pb{E) 



Po 
3 



1^1 



-1/3 



+ \E - fioB\-^/'' + \E + fioBl-^^"") 



(3.10) 



where the three separate contributions come from the 5*^ = 0, ±1 excitations respectively. 
The z = 3 excitations for the 5*^ = ±1 sector are gapped out and this results in a suppression 
of the low energy DOS and consequently also the specific heat and magnetic susceptibilities 



at low temperatures, as shown in Figures 6(a) and 6(b) respectively. As T — )■ 0, only the 
Sz = excitations contribute to the specific heat which thus gets reduced to a third of its zero 
field value. Since no excitations of the z = 3 kind contribute to the magnetic susceptibility 
at T = 0, it is reduced to in comparison to its zero field value. This leads to a suppression 
of the Wilson ratio at low temperatures, as shown in Figure 6(c)[ At temperatures much 
higher than the Zeeman energy, these quantities recover their zero field values. 



c/c, 




T/(/ioB) 




T/(/LJoB) 




T/(/JoB) 



FIG. 6. The suppression of specific heat (a) and the spin susceptibility | (b) | by a magnetic field that 
couples via a Zeeman term and gaps out the Sz = H excitations with z = 3. This also leads to 



the suppression of the Wilson ratio at low temperatures, as shown in (c) At temperatures much 
higher than the Zeeman energy, these quantities recover their zero field values. 



IV. EFFECTS OF WEAK DISORDER 
A. The bond impurity potential 

The Majorana bilinear Hamiltonian does not allow the incorporation of the widely-used 
on-site local impurity model, since 7^ (a;) has to be equal to 1. The simplest kind of spin- 
rotation symmetric local impurity allowed in the spin model is a disrupted bond which 
has the following mean field form, assuming that the disrupted bond is oriented along the 
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direction 5, joining r and r + S: 

Vs{r) = 5 J (5(r) ■ 5(r + S)) -tUj^ b'^b'^ e'^P^"^-"- {e'^-^ - e'"'^) (4.1) 



Here, U (x gSJ is a real number and we have used the (anti) symmetry of the relevant 
operators. 



B. The disorder-averaged self energy in the Born approximation 

The disorder-averaged self energy in the Born approximation is given by (not including 
averaging over different bond directions) 

jj2 f A 1 - cos((p + qr) ■ 6) 

2 [ d^Q 1 ~ cosp • S COS qr ■ <5 + sinp ■ <5 sin qr ■ <5 , , 

= n,^^U J — ^_E^ + ^0+ ^^'^^ 

At this point, we take into account the fact that these bond disruptions are randomly 
oriented in space by averaging the above expression over the three values of 6 related by 
rotations through 27r/3, as shown in Figure 12} We perform this average by using the fact that 



the denominator of the integrand in equation (4.2) is separately invariant under rotations 
of g by 27r/3 and also using the following expressions for averages over 27r/3 rotations over 
the direction of any arbitrary vector d: 

(cosd ■ e)2„,/3 ^ Md) + 0{d') (4.3a) 

(sind ■ e)2„^/3 ^ xe{d) '= cos[3^^] + 0{d') (4.3b) 



In these equations, 9 is the angle between d and e. Using these in equation (4.2), we find 
the rotationally averaged self energy to be'^ 



d^q 1 - (^^(p)(^^(g) + Xx{p)Xx{q) 



27r2 u- Eq + iO+ 

= riimpU^ (/o(w) + fi{uj)^M + f2{uj)Xx{p)) (4.4) 

From this expression we see immediately that as w — )■ 0, both /o,i(w) ~ p{uj) oc uj~^^^ diverge 
as — and thus the Born approximation cannot be justified. This leads us to consider 
the self-consistent Born approximation (SCBA) in the next section. 
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C. The disorder-averaged self energy in the self-consistent Born approximation 
(SCBA) 



The SCBA modifies equation (4.4) to the self-consistent equations 



(c^) - n,^^U y ^-E^- (u;) + zO+ ^^'^^ 

For smaU u and "p these equations can be simplified to: 

4 j 271^ CO -Eg- (w) + Z0+ 

= FoM + Fi(a;y (4.6) 

Using this approximate rotational invariance of Ep*(a;) for small p we can simplify these 
equations to a form that can be easily solved numerically: 

s-,.) . ('i^) r\,, , (4.7) 



4^ / A ./g6_(c^_Ef 



At = 0, these equations can also be analytically solved to the leading order in the disorder 
strength nimpt/^ — > and they yield a purely imaginary value for Y7^^{^). We first deal with 
Fo(0) which depends directly on the momentum cutoff A: 

i^o(O) = f^^) ^dqq't -^^in,^^U') (4.8) 



Att J Jq Att 
The value of -Fi(O) is given by (relabeling Fo(0) by — iFo below, with Tq > 0) 



Numerical solutions to the SCBA equation (4.6) agree with these analytic results (setting 



A ^ 1). The variation of the imaginary parts of the self energy with frequency and impurity 
strength are shown in Figure [7j 
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FIG. 7. Figures showing the variation of the imaginary parts of i^o,i) as defined in equation (4.6) 



and equation (4.7), with frequency as well as the disorder strength. Wherever possible, power law 



fits have been made. In (a) and (b) we have used the disorder strength value Uimp 

[7^ = 5 X 10-^ 



D. The disorder-averaged single particle density of states 



The single particle density of states (DOS) is given by: 



piE) 



Im 



47r3 \E- Ep-J:'p^^{E) 

A~l 

dp 




p 



^p^ - (CO - {Fo{u) + E,{u)p^)y 



(4.10) 



Since Im(Ep*(£' = 0)) is a finite number, we expect the DOS to become constant at low 



energies, instead of diverging as E like in the clean case (3.3). This quenching of the low 



energy divergence in the DOS is shown in Figure 8(a) The variation with disorder strength 



of the low energy saturation value of the DOS is also plotted in |8(b) 
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FIG. 8. Figures showing the variation of the DOS equation (4.10) with energy in |(a)| in the presence 
of disorder with strength nirapU'^ = 5 x 10~^, as well as the variation of the low energy saturation 
value with the disorder strength in |(b)[ 



E. The specific heat in the presence of impurities 



Using the SCBA analysis result that the low energy DOS p(0) is finite, the low temper- 
ature specific heat is found to be 



C 







dT 



6 



(4.11) 



Thus, at very low temperatures, the specific heat is linear in temperature and the coefficient 
of this linear variation is proportional to p(0) and hence to {ni^pU'^)~^^^ (see Figure 8(b)). 



The variation of C/T with temperature is shown in Figure 9(a) where the transition, at 



higher temperatures, to the behavior (3.4) in the clean limit can be seen. 
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FIG. 9. Figures showing the temperature variations of the specific heat divided by temperature 
C/T in |(a)| as well as the static spin susceptibility x in |(b)[ At temperatures larger than a value set 
by the disorder strength nimpU'^ = 5 x 10~^, the quantities regain their behaviors equation (3.4) 



and equation (3.5) in the clean limit. 
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F. The spin susceptibility in the presence of impurities 



The spin susceptibihty is found to be 



Xzz = 2/io 
T^o 2/iop(0) 



dE 



dnF{E) \ 
dE J 
dnF{E)\ /iop(O) 



dE 



(4.12) 



The variation of the spin susceptibihty with temperature is shown in Figure 9(b) where we 
can again see the transition to the behavior (3.5) in the clean limit at higher temperatures. 



G. The Wilson ratio in the presence of impurities — comparison with a 2DEG 



As T — )■ 0, the Wilson ratio for our model is the same as that of a 2DEG of spin 1 
fermions'^ because of the finite DOS at low energies: 



W 



MOP(O) 



Wi 



/2 



7rV(0) 



X 



Tg2DEGT _ 8 _ 1^1 



g2DEG Q 
4 



Wl 



/2 



(4.13) 



where W1/2 is the Wilson ratio of the free spin 1/2 electron gas. We have, as before, assumed 
that the effective magnetic moment of the spins fiQ = fiB, the Bohr magneton. 



Figure 10 shows the variation of this ratio as a function of temperature, showing the 



transition to the clean limit value in equation (3.7) at higher temperatures 



W/Wi/2 

4.2 J — 



10- 



10- 
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FIG. 10. The Wilson ratio W = Tx/C as a multiple of that of the spin 1/2 2DEG W1/2, calculated 
as a function of temperature, in the presence of disorder with strength nimpf/^ = 5 x 10~^. As 
T — )• 0, the value becomes that of a spin 1 2DEG while at larger temperatures it increases to the 



value derived in equation (3.7) for the clean system. 
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H. The thermal conductivity 



The thermal current is^ 



J(g^O,^]^0)= Yl ^''(^ + ^)b-k,-A+n (4.14) 

The thermal conductivity tensor is given bj^ 

^ = - hm hm (4.15) 

where 11 is a tensor whose components are the correlation functions of the thermal current 
components. 11 is diagonal because averaging over three fold rotations makes the off-diagonal 
component J^Jy oc v^Vy vanish. 

For the following calculation, it will be useful to mention these formulae for the quasiparti- 
cle velocities, assuming a low energy long wavelength energy dispersion Eq = tq^ cos{39q + 0) 

cos 9q dq - Eq = 3qh cos(2^q + 0) 

/ cos 6 \ 

Vy = f sin^q dq + de \ Eq = -3qh sin(2^q + 0) (4.16) 

Thus Vq ^ 9g*^, using units in which t = 1. 

The bare thermal polarization bubble (using the renormalized propagators, though) 
yields, after a three- fold rotational averaging that converts v^ y — j- 



p&BZ' ^ ^ 

- /^^/W^°°rfc.(|)'(-^) [(Imact(p,c.))' + (Imact(-p,u;))'] 

y Svv\^^GUP.^ = ^)f (4.17) 



327r3 

T-!>0 9 



327r 

Thus, at low temperatures the thermal conductivity is also linear in temperature (shown in 
Figure 
Figure 



11(b)). The coefficient of this linear variation is plotted vs the disorder strength in 
11(a) We find that the low temperature value of njT varies as (?^impf^^)~^''^ with 
the disorder strength. 

The value of njT in the dmit-131 compound was found to be 0.2 in SI units-, which is 
equivalent to about 332 per triangular spin lattice sheet in units of k\lh^ using the provided 
value of 3nm as the interlayer distance. From Figure 11(a) we find that we require the 
disorder strength to be nimpt/^ ^ 5 x 10^^ to reproduce this value. 
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FIG. 11. The T = thermal coefficient of the thermal conductivity equation (4.17) k is plot' 



ted against the disorder strength in |(a)[ |(b)| shows the variation of this thermal coefficient with 
temperature. 



The vertex corrections to the thermal current can be achieved by replacin^^S] ^j^g quasi- 
particle energy function in the calculation of the quasiparticle velocities by -E'fc+ReSfc(c<;). 
However, ReSfc(a;) is negligible with respect to and so the vertex corrections are negligible 
at low disorder strength^. 



V. CONCLUSIONS 

Let us summarize the characteristic properties of the Majorana spin liquid state on the 
triangular lattice. 

• There are low energy spin excitations near g = which disperse as a; ~ g^. In addition 
there are six Fermi lines which intersect at g = 0, with linear dispersion across the 
Fermi lines. 

• The spin susceptibility, Xi and the specific heat, C are dominated by the z = 2> 
excitations near g = 0; hence x ~ C*/^ ~ T~^/'^, and the Wilson ratio is found to be 
ly -4.2. 

• In the presence of weak disorder, these divergencies saturate at low enough T. Hence 
X ~ C/T ~ T°, and the Wilson ratio W = 8/3, as expected for a Fermi surface of 
S = 1 fermions. 

• The longitudinal thermal conductivity k ~ T as T — with non-zero disorder scat- 
tering. The thermal current is carried mostly by the excitations on the Fermi lines. 

• In the presence of an applied magnetic field, there is no orbital coupling to transverse 
thermal current, and so no thermal Hall effect. 
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• Two-thirds of the lo excitations are gapped out by an apphed magnetic field 

by the Zeeman couphng. Consequently, the specific heat and the spin susceptibility 
are suppressed by the applied field. On the other hand, the thermal conductivity is 
insensitive to the field because it is dominated by the Fermi line excitations, and these 
survive the Zeeman coupling. 

We emphasize that the qualitative aspects of the above results rely only on the assumption 
of a spin liquid ground state on the triangular lattice with SU(2) spin rotation invariance 
and spin-ful Majorana excitations obeying a trivial PSG. 

Many of these properties make our Majorana state an attractive candidate for EtMe3Sb[Pd(dmit)2]2: 
the behavior x ~ C /T ~ k/T ~ T°, and the absence of a thermal Hall effect. An interesting 
distinguishing feature of our theory is that x ^-iid C are suppressed by an applied magnetic 
field, while k/T is not. It would be interesting to test this in future experiments. 
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Appendix A: The relation between the Majorana and spin 1/2 Hilbert spaces 

There are three Majorana operators j^'y^^ per spin. If we consider a collection of an even 
number N of spins, this will lead to a Hilbert space of dimension 2^^/^. We will review here 
how this is equivalent to 2^/^ copies of the 2^-dimensional spin half Hilbert space. 

We assume that the Majorana Hilbert space is composed by randomly picking up pairs 
of Majoranas and representing them by a two-state complex fermion Hilbert space. We can 
show that this space is independent of which scheme of pairing is used. 

Let us now define the following site and bond operators using the Majorana fermion 
operators: 



0{r) = i7^(r)7^(r)7^(r) 



(for ri 7^ r2 only) 



(Ala) 
(Alb) 
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These commute with the spin operators 



[O(ri),S"(r2)] = (A2a) 

0{r) = 27"(r)5"(r) = 25"(r)7"(r) (A2b) 

and also satisfy the following algebraic/ (anti) commutation relations (in the following, 7^ rj 
as long as i ^ j): 

{0(n),0(r,)} = 25,, (A3) 

[r(ri,r2),r(r3,r4)]=0 (A4) 

{r(ri,r),r(r2,r)} = (A5) 

r(ri,r2) = -r(r2,ri); = -1 (A6) 

r(ri,r)r(r,r2)=r(n,r2) (A7) 

[7'^(n), r(r2, ra)] = {r{ri), Tin, r^)} = (A8) 



The consequences of the foregoing relations are as follows. Let us cover the lattice with 
a pattern of bonds (assuming an even number J\f of sites; there are jV!/ (2'^/^(j\/'/2)!) such 
coverings). For such a covering, the J\f/2 number of T operators on the bonds may be 
diagonalized simultaneously and each operator assumes one of the values ±i. There are 
2^/^ such choices. However, for each such choice, the spin operators may be diagonahzed 
simultaneously within each subspace - any spin Hamiltonian is thus not going to mix equiv- 
alent subspaces. Since there were states in the original Majorana Hilbert space, each 
subspace contains 2'^-^/^/2-^/^ — 2^ states which matches the number of states we require 
for the M spin 1/2 's. 

We can explicitly build up such a subspace from the Majorana Hilbert space and demon- 
strate the one-to-one correspondence with the spin states. Let us choose a particular bond 
structure. Now, within every bond, let us label the two sites as m and n. Then, we can 
form the fermion operator 

Cn = ^^^^^f^ (A9) 

which satisfies the usual complex fermion anti-commutation rule (no sum over a) 

Kn. (Cn)^} = 1 (AlO) 



and allows us to define a minimal Hilbert space defined by eigenstates of the number operator 
^mn = {^rnn)'^^mn ~ (l + ^7m7n)/2- the foUowing, wc shall oftcu drop the sub/superscripts 
when there is no ambiguity. 
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The bond operator is proportional to the fermion parity operator specific to that bond 



= i (2n" - 1) {2ny - 1) {2n' - 1) = (AH) 

Using the notation |n^n^n^) = {{c^Y)"' ((c^)^)'^ ((c^)^)" |000), we see that the T = +i 
subspace consists of the states |100),|010), |001) and |111) having odd number of fermions. 
The other case oiT = —i involves states with an even number of fermions. 

In this basis the operators can be expressed as: 

which tells us that the spin 1/2 states are related to the Majorana states in the T = +i 
subspace as follows: 

IS^ ^±1/2,5,;.^ ±1/2) ^W±l^ (^,,,, 
|g; = Tl/2.S,;^±l/2)= (A14b) 

We can proceed similarly and build up a correspondence between the spin 1/2 states and 
a specified subspace, bond-by-bond through the entire collection of spins. We note here that 
the singlet state on a bond is given by the three fermion state |111). 

We can now comment about the relation between the Z2 gauge equivalence apparent in 



equation (1.1) and the exphcit construction of equivalent subspaces above. Such a gauge 
transformation flips the sign of the T operator on the associated bond and thus exchanges 
the subspaces related by flipping that sign. 
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